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Background: A quantitative understanding of neutrino interactions with nuclei is needed for precision era neu¬ 
trino long baseline experiments (MINOS, NOvA, DUNE) which all use nuclear targets. Pion production is the 
dominant reaction channel at the energies of these experiments. 

Purpose: Investigate the influence of nuclear effects on neutrino-induced pion production cross sections and 
compare predictions for pion-production with available data. 

Method: The Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model is used for the description of all incoherent 
channels in neutrino-nucleus reactions. 

Results: Differential cross sections for charged and neutral pion production for the MINERzzA neutrino and 
antineutrino flux are calculated. An estimate for the coherent cross section is obtained from a comparison of data 
with theoretical results for incoherent cross sections. The invariant mass iW) distribution of the A resonances 
produced is analyzed. 

Conclusions: Final state interactions affect the pion kinetic energy spectra significantly. The data for charged 
pion production at MINERi/A are compatible with the results of calculations using elementary data taken from 
an old Argonne National Laboratory experiment. Remaining differences for charged pion production can be 
attributed to coherent production; the data for antineutrino induced neutral pion production, where no coherent 
contribution is present, are reproduced quite well. The analysis of lU-distributions shows that sharp cuts on 
experimentally reconstructed invariant masses lead to shape-distortions of the true W distributions for nuclear 
targets. 


I. INTRODUCTION 

The planned deep underground neutrino experiment 
(DUNE) will operate in a neutrino beam the properties 
of which are similar to the presently used NUMI beam 
at Fermilab. Already ongoing and planned experiments 
in that beam offer the possibility for precision measure¬ 
ments of the differential cross sections for quasielastic 
scattering (QE) and single- and multi-pion production, 
as a function of (reconstructed) neutrino energy. These 
cross sections are essential for an understanding of the 
interactions of the incoming neutrinos with the nuclear 
targets (C, O, Ar) used in these experiments. A quan¬ 
titative understanding of these cross sections is essential 
for the energy reconstruction and - consequently - for a 
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precise extraction of oscillation parameters from observed 
event rates [l| . This extraction has to rely on generators 
and thus the cross sections measured now in experiments 
such as, e.g., MINERzzA could be a basis for a fine-tuning 
of event generators Q for use at the DUNE. 

All events following from a neutrino-nucleus inter¬ 
action can be grouped into two general classes. The 
first one is that of quasielastic scattering (QE), includ¬ 
ing many-body interactions. The second class is con¬ 
nected with pion production, either through resonances, 
t-channel background processes or deep inelastic scatter¬ 
ing (DIS). At the presently running MINER^zA experi¬ 
ment the pion production channels make up for about 2/3 
of the total inclusive cross section while the QE events 
account only for about 1/3 Both event types are 
entangled with each other and cannot be separated by 
purely experimental means. Separation is only possible 
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with the help of an event generator and, therefore, is 
necessarily model-dependent. Recent calculations of the 
inclusive QE response of nuclei 0-0| in their comparison 
with experiment thus have to rely on the accuracy with 
which pion production is described by the generator used 
in extracting the data. 

It is, therefore, essential to obtain a quantitative under¬ 
standing of the dominant reaction channels at MINER^^A 
and DUNE; these are those connected with pion pro¬ 
duction. It is the purpose of the present paper to 
present results of calculations for pion production at the 
MINERVz^A experiment. While first results were al¬ 
ready contained in BQ in the present paper now a 
detailed comparison with MINERi^A data is performed 
using their flux and invariant mass cuts. Over the last 
few years valuable insight into the neutrino-induced pion 
production on nuclei at lower energies has been obtained 
mainly from the MiniBooNE experiment Ill-ll4| . There¬ 


are investigated. 


fore, also the consistency of its data with those obtained 
at the higher energies of the MINERi^A experiment is 
discussed in the present paper. Finally, also some prob¬ 
lems connected with kinematical cuts in the experiment 
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II. METHOD 

The calculations are performed within the transport 
theoretical framework GiBUU [l^ . GiBUU approxi¬ 
mately factorizes the reaction into a very first interaction 
of the incoming neutrino with a bound and Fermi-moviM 
nucleon and the following final state interactions (FSIeJ. 
The latter are described by a numerical implementation 
of the Kadanoff-Baym equations in the gradient approxi¬ 
mation [l3 | , using the Botermans-Malfliet approximation 
for off-shell transport tin. More details about this treat¬ 
ment of final state interactions can be found in Ref. [l^ . 

The groundstate of the target nucleus is assumed to 
be that of a local relativistic Fermi gas with the nucleons 
being bound in a coordinate- and momentum-dependent 
potential that has been fitted to equation of state and 
effective mass data Qi 0 • The hole spectral function 
is given by 


P, 


.{P:E)=g J d^r 0 [pF(r) - |p|] &{E)6 (^E - m*{r, p) -b -b m*^(r, p)^ ; 


( 1 ) 


nucleus 


r 


here PF(r) is the local Fermi momentum given by the 
local Thomas Fermi model and g is a degeneracy fac¬ 
tor. In this spectral function all effects of the nucleon 
potential are assumed to be contained in the effective 
mass m* [0 which depends on location and momentum 
of the nucleon. The corresponding momentum distribu¬ 
tion approximates that obtained in state-of-the-art nu¬ 


clear many-body theory calculations quite well; see Fig. 
4 in 2l| . The initial interaction rates are calculated for a 


nucleon at rest; the results are then boosted to the local 
rest frame of the Fermi-moving target nucleon. 

The single-pion production cross section at fixed neu¬ 
trino energy is then given by (cf. 2^ where all the details 
can be found) 


^ I denote by ’final state interactions’ only the secondary and fol¬ 
lowing collisions. This is different from the nomenclature often 
used in studies of inclusive cross sections for QE where FSI de¬ 


note the potentials felt by the outgoing nucleon in the final state 
of the initial reaction. 
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9 J J [PF(r) - |p|]/corrdCT“®‘^7^PB(r,p)F^(q^,r) . 


( 2 ) 


nucleus 


Here fcorr is a flux correction factor, fcon = {k-p)/{k°p°); 
k and p denote the four-momenta of the neutrino and 
nucleon momentum, respectively. 7^pB(r,p) describes 
the Pauli-blocking and the factor F.n.(q^,r) in ([2]) de¬ 
scribes the effects of all the FSI contained in GiBUU. 
For pions the latter involve elastic and inelastic scat¬ 
tering as well as pion absorption through two-body and 
three-body processes of the types N* + N N + N and 
N* + N + N ^iV-|-iV-|-A^; here N* stands for the A res¬ 
onance and higher excitations of the nucleon. In GiBUU 
these resonances are treated explicitly; they become ex¬ 
cited and then are being propagated until they decay or 
collide with other nucleons; this is described in some de¬ 
tail in [l^. The time-development of the tt — N — A 
dynamics in the nucleus is determined by the resonance 
widths and collision rates alone. In the resonance region 
there is no room for a further free parameter, such as a 
formation time used in other generators 23j. The use 
of formation times, during which interactions of the pro¬ 
duced particle are prohibited, introduces a new parame¬ 
ter and a corresponding arbitrariness into the description 
of pion production. Only in the DIS part, above invariant 
nucleon masses of about 2 GeV, the concept of a forma¬ 
tion time makes sense since it accounts for the widths 
of high lying, no longer separable excitations of the nu¬ 
cleon. Even then, cross sections of the produced particles 
should rise with time until the final hadron has fully been 
formed 241. 


The cross section dcr™®'^ in Eq. stands for the 
pion production inside the nuclear medium. Pions can 
be produced either through nucleon resonances and t- 
channel (background) processes or through DIS. The lat¬ 
ter, denoted in the following by duDiS) is obtained from 
the string-fragmentation model PYTHIA GiBUU 
smoothly switches over been these two pictures around 
an invariant mass of the nucleon of 2 GeV. Below this 


I- 

energy all nucleon resonances are treated with their cor¬ 
rect pion decay branches taken from the PDG 26|. In 
this mass region the pion production cross section con¬ 
sists of a resonance and of a background contribution 
which have to be added coherently in order to obtain 
the full cross section. For the background contribution 
a form taken from an effective field theory treatment of 
pion production up the A resonance region 
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IS 


used. No background contributions for the higher reso¬ 
nances are taken into account. Since transport (or any 
Monte Garlo code) cannot handle the coherence of res¬ 
onance and background amplitudes the coherent sum is 
split up into a resonance cross section and one for the 
background which contains both the squared background 
amplitude and the interference term. Background pions 
are then produced locally, without any time-delay. Be¬ 
cause GiBUU describes only incoherent processes it does 
not contain a coherent pion-production cross section; the 
latter would require a quantum mechanical description. 

The resonance production cross section is given by 

dCT“®‘^ |k'| A“®‘^(p') 


\Mi 


( 3 ) 


dajdll' 327r^ [(fc • p)^ — 

Here M is the nucleon mass, p denotes the nucleon’s four- 
momentum, p' that of the outgoing resonance and k and 
k' that of the initial and final state lepton, respectively. 
The quantities w and fl' give the energy transfer and the 
scattering angle of the outgoing lepton, respectively. The 
in-medium spectral function of the resonance is denoted 
by A“®‘^(p')- 

The cross section for resonance formation m contains 
the square of an invariant matrix element Mji that is ob¬ 
tained by contracting the lepton tensor with the hadron 
tensor. For the latter one has 

= ixr [p + (4) 

where Aajs is (for the A) the spin-3/2 projector and the 
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vertex factor r“^ is given by 

r“^ = — A°‘^\ 7® 


( 5 ) 


for a positive parity resonance. The vector part V is 
taken from the MAID analysis of electron scattering data 
H so that the data for electro excitation of nucleon res- 

I 


= - 


Ct{Q^) 

M2 


onances are reproduced by construction. For the axial 
part A the spin-3/2 transition current contains in prin¬ 
ciple four independent axial form factors but the 
presently available data do not allow to determine them 
separately. Therefore, one of them(C'^) is is set equal to 
zero 




M2 


I 7 


( 6 ) 


with the further simplification = —C^j^ and = 
^M^/((5^ -I- toJ). More details can be found in Ref. 
22 |. 

The absolute strength of the resonance contributions 
determined by C/^(0) is obtained by fitting the available 
pion production data on an elementary target. The two 
datasets available are those obtained at Argonne Na¬ 
tional Lab (ANL) 30| and Brookhaven National Lab 
(BNL) 3l|. These two datasets differ in the relevant en¬ 
ergy regime, with the BNL dataset being higher than the 
ANL one by about 25%. This introduces a correspond¬ 
ing uncertainty into the calculations for nuclear targets. 
Earlier it had been argued that this difference is due to 
a. The author a detailed 

study of the consistency of the various isospin channels 
and the measured da/dQ"^, concluded that the BNL data 
were too high. This has very recently been verified by 
the authors of 33|. In a reanalysis of the old data that 


fixes the flux with the help of the QE cross section it 
was shown that at least for the 7r+ channel the ANL 
data were preferable and also consistent with other data 
from CERN 331. The reanalysis of Wilkinson et al. 331 


seems to settle the question for the correct elementary 
cross section. However, most recently a new theoretical 
calculation of pion production on deuterium has shown 
that even in this small system ESI can play a significant 
role 3^. This then affects the extraction of cross sec¬ 


elementary pion production cross section. 

The extensive pion production data from the Mini- 
BooNE experiment [ll|, [l2| obtained on a CH 2 target 
are consistently higher than the ones calculated within 
GiBUU, both for the ANL and - less so - for the BNL 
input cross sections [l^ . Q . Motivated by the reanalysis 
by Wilkinson et al [33j now the ANL dataset is used as 
default input into GiBUU. In order to illustrate the sensi¬ 
tivity of the results for a nuclear target to the remaining 
uncertainties in the nucleon cross section, however, also 
some results obtained with the BNL input are shown. 

To obtain the pion production cross section the reso¬ 
nance formation cross section is multiplied with the the 
branching ratio for decay into the ttN channel 

dwdD'dH^M dTTdwdD' Ttot 
Here it is assumed that the decay of the resonance hap¬ 
pens isotropically in the rest frame of the resonance. This 
is not a good approximation for an isolated A resonance. 
Therefore, for the ESI process tt-I-A—^A—s-tt-I-N 
a more sophisticated parametrization of the p-wave be¬ 
havior of the A decay is used as explained in detail in 


35j |. Although most of the primary reactions populate 


tions for p and n targets from data obtained with a D 
target. There is thus still some uncertainty left on the 


higher lying states via DIS or higher resonances, the A 
resonance plays quite an essential role in the observable 
pion spectrum as will be discussed later. 

The in-medium effects for pion production are con¬ 
tained both in the spectral function A in Eq. © as well 
as in the branching ratio for the resonance decay into ttN 
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where the final nucleon state may be Pauli-blocked. The 
final state nucleon is bound in a momentum and coor¬ 
dinate dependent potential which - through energy- and 
momentum-dependence - affects the decay width. The 
spectral function contains a collisional broadening which 
is taken from intensive investigations of A properties in¬ 
side the nuclear medium by the Valencia group [^. To 
obtain the correct spectral distribution of pions it is es¬ 
sential to maintain consistency between this collisional 
broadening and the actual collision rates determined by 
the collision cross sections encoded in the generator. 


The theory described so far has been extensively tested 
with the help of photon- 


37l | and electron-induced 


pion production data. This is not the case for the recent 
work of Yu et al. 


which uses the Adler-Nussinov- 


Paschos model for a description of pion production in 
the MiniBooNE and MINER^^A experiments. This model 
does not describe the tt — N — A dynamics correctly and 
indeed the authors of 


themselves ascribe the failure of 


the model to describe the MINERz/A data to the absence 
of all pion-nucleon interactions, except for absorption, in 
that model. 


For the comparison with the experimental neutrino 
data the fixed-energy cross sections described above are 
folded with the flux 


/•lOGeV 

{da)= / dE,<^(E,)(da“'='^(E,.)+daDis(E.)) 

Jl.5GeV 

( 8 ) 

where is the incoming energy distribution (the 

’flux’), normalized to 1. The integration boundaries 1.5 
and 10 GeV are those of the experimental analysis. Fur¬ 
thermore, for the charged-pion cross sections a cut on 
the incoming invariant mass IE < 1.4 GeV has been im¬ 
posed, as in the experimental analysis. All results shown 
in this paper were obtained with the MINERi/A flux for 
a CH target. All cross sections are given per nucleon, i.e. 
the total cross section for GH is divided by 13. 



FIG. 1. (color online) Kinetic energy distributions per nu¬ 
cleon of incoherently produced charged pions on a CH target 
in the MINERi/A neutrino flux. The solid curve gives the 
results obtained with the ANL cross sections as elementary 
input. The dashed-dotted curve gives the same distribution 
calculated with the BNL elementary cross sections. The data 
are from Ref. Ilfll. 



FIG. 2. (color online) Angular distribution of incoherently 
produced charged pions in the MINERz/A experiment (data 
from j^l. The solid curve gives results of a calculation with 
the ANL input, the dash-dotted curve has been obtained with 
the BNL input. 
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III. PION PRODUCTION AT MINER;/A 

A. Charged pion production in the neutrino beam 

In Fig. [T] the calculated kinetic energy distribution for 
the sum of tt'*' and 7r“ is shown in comparison with the 
data [^. The solid curve, calculated with the ANL in¬ 
put, follows the shape reasonably well. The biggest dis¬ 
agreement shows up for kinetic energies < 0.1 GeV; 
for the higher energies the calculation lies at the lower 
end of the error bars so that the disagreement there may 
not be significant. 

Before discussing the difference between theoretical 
and experimental kinetic energy distributions further in 
Fig. [2] the calculated angular distribution is shown in 
comparison with the data. It is evident now that the 
disagreement noticed for the kinetic energy distribution 
is localized at forward angles 9 < 50°. At about 30° the 
calculated cross section amounts to only about 1/2 of the 
measured one while it describes the data quite well for 
9 > 50°. 

The major systematic error in the calculated results 
comes from the uncertainty in the elementary cross sec¬ 
tion. Figures [Hand [2] also show in the dash-dotted curves 
the distributions calculated with the BNL input. As dis¬ 
cussed earlier this is probably an overestimate of the true 
cross section and is shown here only in order to illustrate 
the effects of a change of input. Now the agreement 
between theory and experiment for the kinetic energy 
distribution (Fig. [T]) is better (calculated cross section 
within all error bars). However, the angular distribution 
now becomes considerably worse (see Fig. [5]) with a clear 
overshooting at the intermediate angles around 60°. 

1. Experiment-Theory discrepancies: Coherent Production 

In order to investigate the discrepancies between the¬ 
ory and experiment further the difference between a 
smoothed curve through the data for the kinetic energy 
distribution and the cross section calculated with the 


ANL input is shown in Fig. |3l The error bars shown 
in that figure are only those given by Eberly et al. 
for the data; no systematic error for the GiBUU calcula¬ 
tion has been added. The difference in the kinetic energy 




FIG. 3. (color online)Difference of experimental and calcu¬ 
lated kinetic energy (top) and angular distributions (bottom). 
The error bars are those of the data taken from [^. Note 
that in the upper figure the Tr axis extends down to only 0.04 
GeV. 

distribution shows an interesting structure with a mod¬ 
est peak at about 0.14 GeV kinetic energy and a clear 
peak at the lowest energy of only about 0.04 GeV. The 
cross section amounts to roughly 0.15 ■ 10“^® cm^/GeV 
for energies above about 0.1 GeV; at the energies below 
that the difference rises up to about 0.5-10“®® cm^/GeV 
at 0.04 GeV. Both the modest peak at 0.14 GeV and the 
sharp peak at 0.04 GeV lie above the error bars. 
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For the angular distribution (lower part of Fig. [3]) the 
difference between the calculated and the experimental 
values peaks at about 9 = 18°. At the peak the error 
is about 30%, while it increases to about 70% for angles 
larger than about 9 = 30°. Over all the angular range 
shown here up to about 60° the difference is significantly 
larger than 0. 

GiBUU describes the incoherent pion production quite 
well. This assertion receives some support from the 
fact that it describes the angular distribution of the CC 
antineutrino-induced neutral pion production quite well 
as will be discussed in some more detail later in Sect, 
nr C II The excess at low kinetic energies and small an¬ 
gles, observed here for CC charged pion production, could 
then be due to the coherent production. 

MINERt/A has also measured the coherent cross sec¬ 
tion but without the IFrec < 1-4 CeV cut so 


that these data are not directly comparable with the 
ones extracted here. The measured value without the 
cut amounts to about 3.5 x 10“^®cm^ the integral 


over the difference shown in Fig. [3] amounts to about 
1.9 X 10“^®cm^ and is thus quite reasonable considering 
the different IF-ranges. This is also supported by the ex¬ 
perimental observation that the discrepancy at the lowest 
energy observed between the CENIE prediction and the 
experimental cross section amounts to about a factor of 
1.5 and is nearly entirely due to events with kF > 1.4 
CeV A combined analysis both of the coherent 


and the total 


40(1 data, using the same cuts and generator 


version, would help to clarify the situation. 

Most generators use PCAC based models employing 
the Adler relations The experimental analysis in 


41| finds ’poor agreement’ of the models implemented 


in the standard generators CENIE and NEUT with the 
data; in particular the NEUT version gives a significantly 
broader (and correspondingly) larger coherent cross sec¬ 
tion. Recently, it has been argued that the Adler model 
breaks down for coherent excitations on nuclei ^ both 


at high and at low energies. In a detailed discussion of 
the Adler model it has also been shown that the original 


Rein-Sehgal (RS) as well as the Berger-Sehgal implemen¬ 
tations of that model involve arbitrary assumptions H; 
in addition, the pion-nucleus elastic cross sections are 
not at all described by the RS model (see Eig. 2 in H), 
as implemented in CENIE A complete theory of 


neutrino-induced coherent pion production must involve 
both reaction amplitudes connected with the peripheral 
t-channel scattering on the whole nucleus and, in addi¬ 
tion, those connected with s-channel excitations of the 
nucleon resonances. This is well known from studies of 
coherent electroproduction of neutral pions at higher en¬ 
ergies at JLAB 


47 


where the theory and the data 
exhibit both of these two contributions, one from the 
Primakoff effect and one from individual nucleon inter¬ 
actions. The Adler model belongs to the first category; it 
suffers from the problems just discussed. Models for the 
resonance excitations are still under development 491-151 1. 


2. Comparison with MiniBooNE CC pion production 

The MiniBooNE pion production data required in¬ 
put data even higher than the BNL values 13| so that 
these data are obviously not compatible with the ones 
obtained by MINERz/A. This disagreement between the 
MiniBooNE results and theor y b ecomes larger with in¬ 
creasing for neutrino energies [13|. Sobczyk and Zmuda 
in a detailed comparison of both data sets, have 
recently pointed out that the experimental pion spectra 
of the MINERz/A and the MiniBooNE experiment lie on 
top of each other for kinetic energies above about 125 
MeV. At lower energy the MINERz/A data are consid¬ 
erably higher so that the shapes of the kinetic energy 
distributions in both experiments are quite different. 

For comparison in Fig. |4] also the calculated pion spec¬ 
trum for the MiniBooNE experiment is shown as the 
lowest, dotted curve; for a more detailed presentation 
of CiBUU results for pion production in MiniBooNE see 
Ref. [l3|. This dotted curve should be compared to the 


middle solid line for MINERi/A obtained with a W cut 
at 1.4 CeV. The calculated MINERz/A cross section is 






















roughly by a factor 1.5 larger than that calculated with 
the MiniBooNE flux; this is in agreement with the com- 
jarison of results obtained with the NuWro generator in 
However, the calculated shapes of the spectra are 
quite similar, in contrast to the experimental behavior 
which shows the strong rise at low pion kinetic energies 
in the MiniBooNE experiment. 


Neutrino-induced coherent pion production increases 
with energy and is localized at small kinetic energies. 
It is, therefore, tempting to assign the shape-difference 
to coherent excitation which is expected to be larger 
at MINERj^A. Unfortunately, the theory of neutrino- 
induced coherent pion production is still in a preliminary 
stage (see discussion in the last section). The results 
shown in ^ indicate an increase by about a factor of 2 


between the lower MiniBooNE and the higher MINERVA 
energies. 



FIG. 4. (Kinetic energy spectra for incoherent single charged 
pion production in CC reactions. The solid lines give the re¬ 
sults for neutrino with the ANL input and a cut Wrec <1.4 
GeV. The dashed, uppermost line shows the single-pi pro¬ 
duction spectrum without any tUrec-cut. The dotted, lowest 
curve shows the results for the MiniBooNE flux and a CH 2 
target, again without a Wrec-cut. 



FIG. 5. (color online) Invariant mass distributions of incoher¬ 
ent events in which first a A resonance was excited. The cal¬ 
culations were performed for a C target with the MINERi/A 
flux. Shown are the true VF-distribution without (solid, green 
curve) and with (dashed, red curve) cut. Also shown is the re¬ 
constructed Wiec distribution (red dash-dotted curve) and the 
distributions of the ttN invariant mass Wt,n without (dotted 
green curve) and with cut (red, dot-dot-dashed curve). 


3. Invariant mass cuts 


All the data in Figs. [T] and [2] were obtained with a cut 
on the reconstructed invariant-mass W^ec <1-4 GeV in 
the analysis [^. Here W^ec is defined as 


11/2 ^ -f 2Muj - 


(9) 


where lu is the energy transfer and the four- 
momentum transfer. The experiment does not measure 
these latter two quantities directly but has to reconstruct 
them from a calorimetrically reconstructed incoming neu¬ 
trino energjl^. The cut has been used in to make the 
higher energy data obtained at MINERjzA comparable 
to those obtained earlier at lower energies in the Mini¬ 
BooNE experiment ll| where the population of higher 
invariant masses was energetically suppressed. 


2 It would be interesting to know how much of the final energy 
is actually measured and how much has to be reconstructed by 
using a generator. 
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In order to illustrate the effect of the cut on Wrec Fig.|4] 
also shows the kinetic energy spectrum without any such 
cut (topmost dashed line). This uncut spectrum is larger 
by about a factor of 1.7 at the maximum. The cut on 
Wrec thus has a considerable influence on the published 
experimental cross section. This large effect is worry¬ 
ing since Wrec has no direct physical meaning for bound 
nucleons. 

Wrec is the invariant mass for an interaction with a free, 
unbound nucleon at rest. For a bound and Fermi-moving 
nucleon the correct invariant mass is given by 

W^ = {En + w)2 - (PN + q)2 , (10) 

where Ejy is the energy of the bound nucleon, pn its 
momentum and q the three-momentum transfer. Also 
W is experimentally not directly accessible. 

Both Wrec and W represent entrance channel proper¬ 
ties. For pion production relevant are also the invariant 
mass distributions of the final ttN pairs 

= {En + E^f - {Pt^ + p^f . ( 11 ) 

Here Ej^ and p,r are the energy and the momentum of 
the pion. For values up to about 1.5 GeV W^at contains 
information on the A spectral function and is experimen¬ 
tally directly accessible. 

In order to investigate these different definitions of 
an invariant mass in some more detail various W- 
distributions calculated with GiBUU are shown in Fig. [5] 
There is, first, the true W distribution lEa. [TUll . without 
any cut, shown by a solid green line. This curve has been 
obtained from a calculation without any cut. The same 
true W distribution for a calculation in which the cut 
Wrec <1.4 GeV has been used is shown by the dashed 
red line. This curve shows no sharp cutoff at 1.4 GeV, 
but instead it looses strength on the high W side, start¬ 
ing from the maximunj^ and upwards towards higher W. 


® The peak in Fig. [5] appears at a somewhat lower mass than the 
free A mass; this shift is due to the binding of the A inside the 
nucleus. 


The cut on Wrec, therefore, does not eliminate the pi- 
ons from high-mass excitations, but instead distorts the 
spectral shape of the A. 

The distribution of Wrec (red, dashed-dotted curve) 
looks very different from the true W distributions dis¬ 
cussed so far. It is considerably lower in its peak cross 
section, broader and shifted to higher masses until it is 
cut off at 1.4 GeV. It shows no resemblance to the A 
spectral function. 

Finally shown are also, by the green dotted and the 
red dot-dot-dashed lines, the invariant mass distributions 
of the final nN pairs. Again, the red dot-dot-dashed 
curve obtained in a calculation with a cut on Wrec misses 
strength from 1.15 GeV on upwards. No sharp cut in 
WttN appears. Since the events in both distributions are 
subject to strong final state interactions they are lower 
than all the others. 

B. Charged pion production in the antineutrino 

beam 

For completeness in Fig. |6]the charged pion spectrum 
(sum of 7r“ and 7r+) in an antineutrino beam is shown. 
These results were obtained with a flux 1.5 GeV < E^ < 
10 GeV and no invariant mass cut. The shape is very 
similar to that obtained for the neutrino beam (topmost 
curve in Fig. II, but is lower by about a factor of 2 due 
to the different V — A interference for antineutrinos. In 
Fig. [7] the angular distribution is given for the incoherent 
single pion production in the MINERVi/A antineutrino 
beam. 

C. Neutral pion production in the neutrino and 

antineutrino beams 

The 7r° production both for neutrino and antineu¬ 
trino induced CC reactions is interesting because there 
is no coherent excitation in this channel so that these re¬ 
sults should be directly comparable to experiment. The 
GiBUU results discussed in this section also obtained 
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FIG. 6. Kinetic energy spectra for incoherent single charged 
pion production in CC reactions with the MINERVi/A an¬ 
tineutrino beam between 1.5 and 10 GeV. No VFrec-cut has 
been used. 



FIG. 7. Angular distribution of incoherently produced 
charged pions in the MINER/^A experiment with an antineu¬ 
trino beam between 1.5 and 10 GeV. No IFrec-cut has been 
used. 

with a flux 1.5 GeV< < 10 GeV and no invariant 
mass cut using the ANL input. 

The two upper curves in Fig. [8] give the results for 
neutrino-induced and the two lower ones for antineutrino- 
induced neutral pion production. The effects of FSI are 
remarkable: the cross section after FSI are significantly 
larger than those before, in particular around 0.1 GeV 
kinetic energy. Furthermore the shape is significantly 
distorted by FSI. The strong overshoot at 0.08 GeV is 


FIG. 8. (color online) Kinetic energy spectra of incoherently 
produced single for CG reactions on a CH target using 
the MINERVi^A neutrino and antineutrino fluxes between 1.5 
and 10 GeV. The solid lines give the results for neutrino (up¬ 
per, green) and for antineutrino (lower, red) beams with FSI 
included. The dashed lines, labeled with ’inch, give the cor¬ 
responding results before FSI. No Wtbc cut has been used. 



FIG. 9. (color online) Angular distribution of incoherently 
produced single vr® in the MINERz/A experiment with the 
neutrino (green upper line) and antineutrino (red lower line) 
fluxes between 1.5 and 10 GeV. No Wrec cut has been used. 

followed by an undershoot at around 0.24 GeV. For some¬ 
what higher kinetic energies the cross section is again 
increased by the FSI and only from about 1 GeV on 
upwards there is a slight attenuation by final state in¬ 
teractions. The cross section after FSI shows a net in¬ 
crease compared to that before FSI. This increase is due 
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to charge-transfer reactions from the dominant chan¬ 
nel in the i^-induced production and the 7r“ channel in 
the V induced one. The importance of this charge trans¬ 
fer was already pointed out in [ 23 ; here it should be 
noted that GiBUU has been extensively tested against 
pion charge-exchange reactions on nuclei 35, H- 


The region around T^r = 0.24 GeV reflects the strong 
pion absorption through the pion-less decay of the A res¬ 
onance. This indicates that even at the high energies of 
the incoming neutrino beam in the GeV region, the A 
resonance plays an essential role in the FSI. Initially pro¬ 
duced energetical pions cascade down through a sequence 
of elastic or inelastic scattering, possibly connected with 
charge transfer. The pions that finally end up in the A 
region can then be absorbed. Those pions that are slowed 
down even further can no longer be absorbed into a A 
thus causing the strong peak at 0.08 GeV. 

The angular distributions for 7r° production shown in 
Fig. 0 look very similar to those for charged pion produc¬ 
tion. 


1. Neutral pion production at MINERvA 


The calculated results for CC antineutrino-induced 
neutral pion production can be directly compared with 
data 53 that became available about 1 week after af¬ 
ter a first version of thepresent paper paper had been 
uploaded to the arXiv The main difference to the 


charged pion data, besides the different flavor of the in¬ 
coming beam, consists in selecting the incoming neutrino 
energy to be between 1.5 and 20 GeV (instead of 10 GeV 
for the charged pions) and the absence of any W cut. 
The absence of a IV-cut has a more significant effect on 
the cross sections than the raising of the upper neutrino 
energy; the latter has about a 10% effect. The results 
obtained with these these experimental specifications are 
shown in Figs. ITUlandfTTl Because of the absence of any 
GG coherent 7r° production of neutral pions these cal¬ 
culations are directly comparable with experiment. The 
comparison of the experimental and theoretical angular 



FIG. 10. (color online) Kinetic energy spectra of incoherently 
produced single for CC reactions on a CH target using the 
MINERVi^A antineutrino flux between 1.5 and 20 GeV. Data 
are from [^, converted into da/dT-,^. 



FIG. 11. (color online) Angular distribution of incoherently 
produced single in the MINERi/A experiment with the 
antineutrino flux between 1.5 and 20 GeV. Data are from 



distribution in Fig. 1111 shows a very good agreement be¬ 
tween both. In particular, there is no excess at small 
angles, as it was the case for the charged pion produc¬ 
tion discussed in earlier sections. This reflects the ab¬ 
sence of any coherent production. The comparison of 
the calculated kinetic energy distribution with the data 
still exhibits an underprediction at the lowest kinetic en¬ 
ergies around 40 and 80 MeV. On the other hand, for 
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kinetic energies above abont 0.3 GeV the data are lower 
than the theoretical prediction. 

This behavior lends some support to the earlier discus¬ 
sions of a coherent contribution to the measured charged 
pion cross sections in Sect. IIII A II The overshoot of the 
kinetic energy distribution at higher energies, combined 
with the underprediction at the lowest energies, may in¬ 
dicate an underestimate of FSI charge exchange reactions 
for the pion and/or problems in the experimental back¬ 
ground subtraction. 


IV. SUMMARY AND CONCLUSION 


In the energy regime of the planned DUNE pion pro¬ 
duction, either through resonances or DIS, represents the 
dominant reaction channel. It is, therefore, mandatory 
to obtain a quantitative understanding of this process. 
Quasi-elastic scattering, which has received a lot of the¬ 
oretical attention recently, is necessarily entangled with 
pion production in the actual observables and can not be 
separated without the help of an event generator. Thus, 
the quality of the QE data can never be better than the 
quality of the generator used to remove the pion contam¬ 
ination. On the other hand, pion production initiated by 
neutrino interactions is interesting in itself since it may 
provide information on the axial couplings of resonances. 

In this paper primarily the recent MINERj^A data on 
charged pion production were analyzed. Unfortu¬ 
nately, the crucial input to such calculations, the pion 
production cross section on isolated nucleons, is still 
somewhat uncertain. While the discrepancy between two 
older data sets, from ANL and BNL, has recently found 
a reasonable explanation in terms of flux uncertainties in 
the BNL experiment 33|, there are still lingering prob¬ 


lems with the extraction of nucleon cross sections from 


experiments using deuterium targets [3^. This latter 
problem clearly deserves more theoretical work, but ulti¬ 
mately a dedicated experiment using hydrogen targets is 
needed to clarify this point. 

The calculations for charged pion production repro¬ 


duce the measured kinetic energy spectra quite well, ex¬ 
cept for the lowest kinetic energies around 100 MeV, 
where the calculation comes out significantly lower than 
the experiment. This discrepancy can be localized in 
events connected with small scattering angles. At for¬ 
ward angles, the calculation significantly underestimates 
the experimental cross sections. These discrepancies, 
both in kinetic energy and angle, are compatible with 
contributions expected from a coherent excitation pro¬ 
cess. A reanalysis of the experimental MINER^^A data 
on coherent pion production , using exactly the same 


0 , 


cuts and generator tune as the Eberly et al analysis 
would help to verify this explanation. This explanation 
also finds some support in the observation that the angu¬ 
lar distribution of the charged current neutral pion pro¬ 
duction data in the antineutrino beam 5^, where no 


coherent component can be present, is reproduced quite 
well by the present calculations. 

The comparison of the MiniBooNE pion production 
data with those obtained by MINERzzA comes to the 
same conclusion as the recent analysis in 5^ that these 
two datasets are not compatible. The shape discrepancy 
may, however, reflect the presence of a coherent compo¬ 
nent in the spectra which will be larger at the higher 
energy of the MINERi/A experiment. It is hoped that 
pion production data from the T2K experiment can shed 


0 . 


some light on this problem 

In an earlier paper [0 the authors had already noted 
that the flux cut used in the MINERz/A analyses intro¬ 
duces a model dependence into the comparison of theory 
with data. The cuts used in theoretical calculations are 
those for true energies whereas the experimental cuts can 
only be done for reconstructed energies. It was shown in 
[lOj that the flux cuts used by the MINERz/A experi¬ 
ment can significantly distort the event rates, requiring a 
fairly large and generator-dependent correction through 
energy-migration matrices. Here now the effects of cuts 
imposed also on the reconstructed invariant mass tUrec 
in obtaining the pion production data were investigated. 
The cut has a major influence on the total pion produc- 
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tion cross section. This was to be expected since a large 
part of the pions is produced by DIS events. A cut at 
1.4 GeV, as employed in the experiment, was intended to 
remove these contributions and enrich the A resonance 
contribution. This, however, is not what the Wrec cut 
actually does. Instead it cuts off strength over a large, 
high-mass part of the A spectral function, starting al¬ 
ready around its peak value. The explanation for this 
behavior lies in the Fermi-motion and binding energy of 
the nucleons inside the nuclear target. A comparison of 
pion cross sections obtained with this cutoff on Wrec with 
data in other experiments and at lower energies is then 
difficult. In addition, this cut - being an entrance channel 
cut - has to rely on the reconstruction of neutrino energy 
which introduces a model dependence into the data. A 
cut on the n — N invariant mass in the final state would 
be free of this problem and is thus preferable. 
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